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Physical Effects in Cavitating Flows 
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A discussion of the physical effects in cavitation naturally divides 
itself into a part concerned with microscopic aspects and a part concerned 
with macroscopic aspects. The microscopic features of cavitating flows 
are related to the properties of nuclei in liquids and to the moderate 
tensile strengths which are usually encountered in flows. The macroscopic 
features are concerned with the growth of vapor or gaseous cavities from a 
small initial size and with their eventual collapse. The growth and collapse 
of vapor or gaseous cavities are macroscopic problems of a transient 
nature. There are macroscopic problems of a quasi-steady nature which 
are encountered in cavity flows — or in supercavitating flows, as they are 
sometimes called. We shall not consider these quasi-steady flows in any 
detail here. 

TENSILE STRENGTH OF LIQUIDS 

As is well known, the theoretical tensile strength of a liquid free of 
contaminants is very large — much larger than is ever observed. The 
theory of this tensile strength is statistical in nature so that one deter- 
mines the probability of appearance of a small vapor bubble under a given 
tension (ref. 1). The formation of a, bubble of radius Ro requires an in- 
crease in free energy of 47rffo 2 <r where <j is the surface tension. As was 
recognized a long time ago by Gibbs, the pressure within the bubble 
exceeds the pressure in the liquid by 2<r / Ro so that the net increment in 
free energy is reduced to 4ir/f 0 2 <j/3. The probability, P, of appearance of a 
bubble is dominated by this free energy : 

P^g-irR^/mT) ( 1 ) 

In equation (1), k is the Boltzmann constant and T is the absolute 
temperature. It is apparent that this probability is very small even for 
large tensions. For example, with <r = 72 dynes/cm (water at 27° C), 
lnP=— 73 when R tJ = 10~ 7 cm; the corresponding tension is 1400 atm. 
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While there is some reduction in the macroscopic value of the surface 
tension constant when the radius of curvature becomes very small (ref. 2), 
this reduction is not significant unless Ri, is less than 10 -7 cm. We must 
conclude, therefore, that a homogeneous, pure liquid should have tensile 
strengths of the order of several thousand atmospheres. It is pertinent to 
remark that this large tensile strength would be predicted for liquids con- 
taining dissolved gases, so long, of course, as the gas is truly dissolved — 
that is to say, so long as the gas is distributed through the liquid on a 
molecular scale. 

The conventional resolution of the disagreement between the theory of 
the tensile strength of a liquid and the experimental findings of much 
lower values is made with the introduction of “nuclei” in the liquid. 
Various properties of these nuclei have been proposed and we shall wish 
to discuss some of these here, but their function is to provide sites at 
which bubbles may grow from some effective size, Ro, which is large 
enough to give a moderate value of 2cr /Ro. In the statistical theory of the 
strength of a pure liquid, a cavity of radius Ro must be formed where no 
nucleus exists through a statistical fluctuation. Such fluctuations are 
most improbable unless the spherical cavity is of the order of molecular 
dimensions. If the cavity in the liquid already exists at a larger radius, the 
corresponding tensile strength will be reduced. 

One model which has been suggested for nuclei is that they are solid 
particles such as dust particles. It is easy to show (ref. 1) that such solid 
particles must not be wetted by the liquid if they are to serve as sites at 
which a macroscopic cavity is to grow. For water one would expect an 
abundance of hydrophobic dust particles, and, to give agreement with 
observations, “ordinary” water should have such particles with sizes at 
least as large as /£o~10~ 4 cm; very carefully prepared “clean” water 
should have particles no larger than f? 0 ~10~ 5 cm or 10 -6 cm. 

Solid particles will tend to settle out from the liquid so that there should 
be some tendency for the tensile strength of ordinary water, for example, 
to rise with time if the sample is quiescent and protected from the at- 
mosphere. Such an “aging” effect has been observed (ref. 3). It must be 
admitted, however, that this aging effect could be explained by un- 
dissolved air bubbles rising out of the liquid. The tensile strengths of very 
clean liquids should not show a significant aging effect since very small 
nuclei, say i2o~10 -6 cm, can be maintained in suspension by the mecha- 
nism of Brownian motion. In a macroscopic flow, much larger particles 
will be effectively maintained in the liquid since their settling velocities 
are so small. 

So far we have not considered any effects of dissolved gases on the 
tensile strength of the liquid. In the application of a single tension pulse, 
or of a steady tension, this picture of the tensile strength as being deter- 
mined by suspended, solid nuclei would indicate that there should be no 


PHYSICAL EFFECTS IN CAVITATING FLOWS 343 


significant effect from dissolved gases. An increase in the effective size of 
solid nuclei by layers of gas adsorbed from solution would be unimportant 
since such a layer would have a thickness determined by the range of 
intermolecular forces. This range is of the order of 10 -7 cm. There is, 
however, a macroscopic mechanism whereby dissolved gases can play a 
significant role; this mechanism is the surface tension resultant at the 
interface of three media: liquid, gas, and solid. The surface tension forces 
can be annulled, or even reversed, for such combinations so that un- 
dissolved pockets of gas can be stabilized in crevices of solid nuclei. 

Undissolved bubbles of gas have been considered as possible nuclei 
from which macroscopic cavities can grow under reduced pressure. Such 
gas bubbles are, however, unstable. In an undersaturated solution the 
bubbles dissolve; even in a saturated solution gas bubbles will dissolve 
because of surface tension (ref. 4). In a sufficiently supersaturated solu- 
tion, gas bubbles will grow indefinitely. While it must be acknowledged 
that the mass diffusion process which leads to collapse or growth of these 
gas bubbles is very slow compared with the dynamic process of vapor- 
bubble growth or collapse, it is difficult to explain a long-term modest 
tensile strength by means of undissolved gas bubbles unless they are 
attached to a solid surface. Such attachment of air bubbles with dimen- 
sions less than 10~ 4 cm has indeed been observed to exist on a rather 
permanent basis on hydrophobic particles in water (ref. 5) . 

There is a mechanism whereby the time required for the dissolution of 
very small free bubbles might be increased. Impurities, dissolved or dis- 
tributed in the liquid, could tend to deposit on the surfaces of gas bubbles. 
This process would decrease the effective surface tension constant of the 
gas/liquid interface and would thereby represent a configuration of greater 
stability since the free energy would be decreased. The decrease in the 
surface tension w r ould in itself lengthen the time for disappearance of the 
bubble and, in addition, could decrease the diffusion rate of gas through 
the bubble boundary into the liquid. An anomalous behavior in small 
oscillating gas bubbles has been observed (ref. 6) which might be due to 
such an adsorbed film. Some anomalous persistence of small bubbles has 
also been claimed (ref. 7) and might be similarly explained. 

A liquid should show the same tensile strength for a single pulse of 
tension or for any nonrepetitive tension regardless of whether the nucleus 
is a solid, nonwetted particle or happens to be a small gas bubble. The 
growth from a nucleus to a macroscopic cavity under these conditions is 
so rapid that diffusion of dissolved permanent gas into the bubble from 
the solution does not take place. 

If there are sufficiently long times available, diffusion of dissolved gases 
can play a central role in increasing the size of nuclei. Such long times are 
produced when a liquid is exposed to oscillating pressure fields. If the 
liquid has a concentration of dissolved gas, c w which is the equilibrium 
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value for the liquid pressure P 0 , then a bubble of radius R 0 can gain mass 
when the liquid pressure undergoes periodic oscillations with angular 
frequency co: 

P(<) =P 0 (l+« sin id) 

This process is known as rectified diffusion, and the rate of increase of 
mass in the bubble averaged over many oscillation cycles is (ref. 8) 


dm 

dt 


8ir 

= — Dc K R 0 f 

O 


where D is the coefficient of diffusivity of the gas in the liquid. Clearly, 
rectified diffusion can stabilize a gas bubble which would otherwise dis- 
solve. A gas bubble of radius R 0 in a liquid at pressure Po, with equilib- 
rium-dissolved concentration c w will lose mass at the following rate 
(ref. 4) : 




The threshold pressure amplitude for rectified diffusion which would just 
balance this loss is 


AP = P max — Po = 



If A P exceeds this value, the bubble will grow — but at the slow rate 
characteristic of diffusion processes. 

Another situation in which the long times required for the diffusion 
process are available is familiar in the flow of water over a submerged 
body so that the static pressure is reduced. The liquid may then be 
appreciably supersaturated, and a small bubble may grow on the surface. 
Such bubble growth has been observed in the boundary layer of a model 
in the working section of a water tunnel. When the bubble has become 
large enough, it leaves the surface of the body, enters the flow, and may 
serve as the nucleus for formation of a macroscopic vapor bubble. 

It may be of interest to note that Einstein gave a theoretical connection 
between gas diffusivity in a liquid and the liquid viscosity : 


(jirfia 

where a is the radius of the diffusing molecule, T is the absolute tempera- 
ture, and fi is the coefficient of viscosity. From this relation it is evident 
that a highly viscous liquid like glycerol should have a coefficient of 
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diffusion which is approximately 10~ s times that of water. The implica- 
tions of such a small diffusion coefficient in a viscous liquid for gaseous 
cavitation are indeed interesting. 

To summarize our view of the effective tensile strength of liquids, the 
small observed tensile strengths are due to the presence of nuclei which 
should for the most part, consist of nonwetted solid particles. Dissolved 
gases should play only a secondary part in tensile strength and cavitation 
except where slow processes can be effective, and typical of these is 
rectified diffusion with oscillating pressures or the flow of supersaturated 
dissolved gas into a bubble attached to a solid. 


DYNAMICS OF CAVITATION BUBBLES 


Growth of Vapor Bubbles 

In the growth of a vapor bubble, we may suppose that the rate of 
growth gives bubble boundary velocities which are small compared with 
the speed of sound in the vapor. It follows that the pressure within the 
bubble will be close to the equilibrium vapor pressure corresponding to 
the temperature in the liquid at the boundary. As the bubble grows, the 
latent heat of evaporation must be supplied at the liquid boundary. The 
associated cooling will reduce the vapor pressure and may therefore 
affect the rate of growth. We may easily find the conditions under which 
this cooling effect becomes important. 

Let us suppose that a vapor bubble grows to a macroscopic size, R, in a 
time r from some initial size R<><KR. Then the total mass of vapor which 
is evaporated into the bubble is (47r/3).fiV where p is the vapor density. 
If L is the latent heat of evaporation per gram, then the total heat re- 
quired is 

Q= (4*-/3)fiyL 

This heat is taken out of a liquid layer around the bubble which has an 
effective thickness determined by the thermal diffusivity of the liquid, k, 

k 

K = — 

pc 

where k is the thermal conductivity, p is the density, and c is the specific 
heat of the liquid. The thickness of the liquid layer is given approxi- 
mately by 

d^((cr ) 1/2 

We see that the estimate of the thickness is not sensitive to the choice 
of r. The volume of the liquid layer is then of the order of magnitude 
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4jt R 2 d, and the corresponding mass of liquid is 4irfl 2 p d. The temperature 
drop in this mass of liquid, AT, is then 

m (4 t/3 )R*p'L R p'L 
4ir R 2 pcd 3 d p c 

The drop in vapor pressure, A p v , associated with this temperature de- 
crease may be estimated from the Clausius-Clapeyron equation, which 
gives 

dp t LM 
dT BT* 


where M is the molecular weight of the vapor, B is the perfect gas con- 
stant, and T is the absolute temperature. This expression for the rate of 
change of vapor pressure with temperature comes from supposing that 
the vapor behaves like an ideal gas and, further, that p5>p ' . Both these 
approximations are fairly accurate. It follows that 



AT 


and finally 


A Pv 

Pv 


LM R p'L 


( 2 ) 


BT 2 3d p c 

The condition for a negligible effect from cooling by evaporation is that 

^«i 


When this condition is satisfied, thermodynamic effects are unimportant 
in cavitation bubble growth. We then say that we have bubble growth 
in a “cold” liquid since it is evident that both p' and ( dp v /dT ) decrease 
rapidly with temperature. 

In a cold liquid, since thermodynamic effects are unimportant, we may 
disregard the energy equation and consider only the momentum equation. 
Since we have already supposed that the bubble growth rate is small 
compared with the speed of sound in the vapor, it will certainly be slow 
compared with the speed of sound in the liquid. The flow in the liquid is 
then essentially incompressible, and if, in addition, we take the growth to 
be spherically symmetric, the momentum equation integrates to the 
Bernoulli equation. From the Bernoulli equation one gets the equation of 
motion of the bubble boundary 


p(R) —Po 
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where P 0 is the pressure in the liquid at infinity and 'p(R) is the pressure 
in the liquid at the bubble boundary. The boundary condition of con- 
tinuity of stress across the bubble boundary gives 

2(7 it 

- 4/1- (4) 

The problem of spherical growth in a cold liquid is thus completely deter- 
mined. We may observe that the left side of equation (3) may be written 
in the form 


1 d 
2 R 2 R dt 


(R 3 R 2 ) 


and one can easily get the asymptotic growth velocity for fixed P<>: 

Hi^r 

since for large R the terms 2 cr/R and 4 fiR/R in p (R) become unimportant. 

When the condition A is not met, the dynamics of vapor-bubble 
growth are strongly coupled to the energy equation, or the heat-flow 
equation. The physical situation now corresponds to tension in a “hot” 
liquid. We consider, for simplicity, the case in which the static pressure at 
infinity, P 0 , is constant, and we assume further that the temperature in 
the liquid is To. The liquid must be superheated if a vapor bubble is to 
grow; that is, the vapor pressure at the liquid temperature, p.(T 0 ), 
exceeds the ambient pressure. We introduce the temperature Tb, which is 
the “boiling” temperature for the ambient pressure P 0 ; that is, Tb is the 
temperature such that p v (T b ) — Pa- Then 


p„(To) >p v (Tb) =Po 

We characterize the “hot” liquid case as that in which the cooling effect 
lowers the temperature in the bubble so that it approaches Tb as the 
bubble grows. When the cooling effect is so large, we have the limiting 
case in which inertial effects are unimportant so that the growth is con- 
trolled entirely by the rate of heat flow. Clearly, the cooling effect can 
never lower the temperature, T, within the bubble so far that p„(T) <P 0 . 
Growth cannot then proceed, and the heat inflow would cease together 
with any cooling. In the asymptotic case, we can expect that the tem- 
perature T within the bubble would exceed Tb by only a small amount, 
so that we may estimate the heat inflow requirement approximately by 
taking Tc^Tb- The temperature in the liquid rises from its value at the 
bubble wall to T 0 in a distance of the order of magnitude (d) 1/2 at the 
time t. This distance is the heat diffusion length. The temperature 
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gradient in the liquid, then, is approximately (T 0 —T b )/(id) 1/2 , and the 
heat inflow into the bubble is 




(T 0 -T b ) 

(*0 I/2 


( 6 ) 


where R = R(t). Now the heat requirement per unit time for evaporation is 


Q=L 


L 

dt 


(frcR*p')~4*R 2 RLp' 


( 7 ) 


If we equate these two values of Q, we find the asymptotic growth speed 


* k ( T 0 -T b ) 

ff/N/ ■■ 

~Lp' (* 0 1/2 


( 8 ) 


We see that the two limiting cases for growing vapor bubbles have very 
simple behavior described by equations (5) and (8). The intermediate 
situation of a “lukewarm” liquid does not have this simplicity, although 
we can be assured that the R (<) trajectory must he between those specified 
by the limiting curves just described. An accurate solution of the general 
case has analytical complications because the heat-flow equation is 
coupled to the momentum equation. The essential complication comes 
from the heat flow across a moving boundary R(t), which is implicit in 
the problem (ref. 9) . 


Collapse of Vapor Bubbles 

When the ambient pressure Pa exceeds the vapor pressure p„, so long as 
these pressures remain constant we may use the momentum equation (3) 
to determine the collapse velocity. If we suppose, for the moment, that p v 
remains constant, and if, in addition, we neglect the effect of viscosity, 
then we easily find that 


R 2 = 


2(Pp — p„) 
3p 


/w v WPo’X 

\R* ) P R \R 2 / 


It follows that R increases in magnitude like Rr 3/2 as R approaches zero. 
Clearly, many of the simplifying assumptions which have been made 
cannot remain valid toward the end of the collapse motion. 

Even in a cold liquid, a stage in the collapse will be reached in which the 
heat of condensation will produce a temperature rise so that p v will in- 
crease. In this case, however, it may be shown by numerical treatment 
that the time for collapse and the trajectory for R(t) are not greatly 
affected by this temperature rise. An additional effect which is of particu- 
lar significance for cavitation or cold bubble collapse is the effect of 
compressibility, which has been studied by many authors. 
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There are evidently many complications which appear in the theory of 
the collapse in addition to the effect of compressibility. Among these are 
the instability in the spherical shape and the difficulty of determining the 
equation of state of the compressed vapor. All these difficulties appear 
only in the last stages of the collapse, when the radius is quite small. It is 
known, however, that a physical ease of particular interest is that of the 
collapse of cavitation bubbles near a solid boundary. For this physical 
situation, the presence of a neighboring solid boundary changes the 
collapse behavior early in the bubble history before the difficulties char- 
acteristic of the late stage of collapse can appear. The bubble is deformed, 
and a jet is formed directed toward the solid boundary (ref. 10). This 
result not only simplifies the theoretical analysis by making it possible to 
avoid these difficulties, but also has important implications for the 
mechanism of cavitation damage. 


REFERENCES 

1. Plesset, M. S., Tensile Strength of Liquids. Cavitation State of Knowledge, ASME 

Symposium, June 1969, p. 15. 

2. Kirkwood, J. G., and F. P. Buff. Statistical Mechanical Theory of Surface 

Tension. J. Chem. Phys., Vol. 17, 1949, p. 338. 

3. Strasbero, M., Onset of Ultrasonic Cavitation in Tap Water. J. Acoust. Soc. 

Am., Vol. 31, 1959, p. 163. 

4. Epstein, P. S., and M. S. Plesset, On the Stability of Gas Bubbles in Liquid- 

Gas Solutions. J. Chem. Phys., Vol. 18, 1950, p. 1505. 

5. Liebermann, L., Air Bubbles in Water. J. Appl. Phys., Vol. 28, 1957, p. 205. 

6. Exner, M. S., and W. Hampe, Air Bubbles in Water. Acoustica, Vol. 3, 1953, 

p. 17. 

7. Turner, W. R., Microbubble Persistence in Fresh Water. J. Acoust. Soc. Am., 

Vol. 33, 1961, p. 1223. 

8. Hsieh, D-Y., and M. S. Plesset, Theory of Rectified Diffusion of Mass into Gas 

Bubbles. J. Acoust. Soc. Am., Vol. 33, 1961, p. 206. 

9. Plesset, M. S., and S. A. Zwick, The Growth of Vapor Bubbles in Superheated 

Liquids. J. Appl. Phys., Vol. 25, 1954, p. 493. 

10. Plesset, M. S., and R. B. Chapman, Collapse of a Spherical Vapor Cavity Near 
a Solid Boundary. J. Fluid Mech., Vol. 47, 1971, p. 283. 



350 CAVITATION 


DISCUSSION 


J. M. ROBERTSON (University of Illinois) : In turbomachinery flows 
with liquids, cavitation may enter to cause noise, loss of machine effective- 
ness, vibration, or material erosion damage. The designer is faced with 
the questions of when and where cavitation will appear, how the flow is 
modified, and the nature of the cavity collapse. A clear and lucid picture 
has been given of the matter of tensile strength of liquids relevant to 
cavitation occurrence and to subsequent consideration of the transient 
dynamics of cavitation bubbles, both growth and collapse, as related to 
damage. And, through concentration on the basic physics, this presenta- 
tion has been achieved without recourse to a single figure or complex 
mathematics. 

In the matter of tensile strength, the author’s review brings this matter 
down from the large values of pure liquids to the rather small values 
usually encountered in engineering. Concerning “dissolved gas” effects, 
I request the author’s interpretation of the rather well-known experiments 
in which pressurization of nonpurified water led to much larger tensile 
strengths. My rough interpretation is that although the gas in the sample 
initially was not visually evident, it was only dissolved on pressurization. 

The author suggests that often for cavitation to appear with oscillation, 
and in body flows, time must be available for gas diffusion into bubbles 
to occur. He notes that with highly viscous glycerol the diffusion would be 
significantly less than with water, thus implying a delay in cavitation 
occurrence. Cavitation experiments have dealt with glycerol, but concern 
was not with the initial appearance of cavitation. An even more viscous 
fluid (pluracol) was used by S. Sundaram (ref. D-l) for studies of 
cavitation inception (actually desinence) for a friction-type valve. The 
cavitation number was found to decrease by a factor of two with a 
doubling of the Prandtl number, which is directly proportional to 
viscosity. 

In view of the use of fluids other than water in some machines, the 
discussion of bubble growth in “hot” liquids versus that in “cold” liquids 
and the resulting thermodynamic effect is germane to a host of turbo- 
machinery flows. In an extensive review of cavitation damage, A. 
Thiruvengadam (ref. D-2) made a suggestive comparison of several 
fluids on this basis. 

R. HICKLING (General Motors Research Laboratories) : Professor 
Plesset has given a succinct summary of our present understanding of 
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some of the physical effects involved in the formation, growth, and 
collapse of bubbles. The theory that he presents allows a clear con- 
ceptualization of the processes involved. 

The assumption of a uniform equilibrium behavior of the vapor inside 
the bubble permits considerable simplifications of the theory, but, in some 
cases, these may not be entirely realistic. It should be emphasized here 
that by a nonequilibrium behavior we mean deviations from the Clausius- 
Clapeyron equation — i.e., the assumption that the liquid and vapor 
phases are continuously in equilibrium with each other. Calculations for a 
bubble containing a permanent gas (ref. D— 3) have shown that significant 
thermal gradients develop in the gas next to the cavity wall, even when 
the bubble wall is moving relatively slowly. A similar situation is expected 
to occur with a vapor. Such nonuniformity of temperature naturally has 
associated with it a nonequilibrium behavior in the vapor, and, under 
certain conditions, the interaction could become quite pronounced. 

It has been argued that, because the thermal diffusivity of the vapor 
in the bubble is much greater than that of the surrounding liquid, the 
liquid presents a kind of thermal barrier to the flow of heat in and out of 
the bubble. It is reasoned that the vapor responds relatively quickly to 
thermal changes and consequently has a more uniform temperature dis- 
tribution, while the liquid responds relatively slowly and has steep 
temperature gradients in the region close to the bubble wall. Surely this 
argument does not completely prove its case. It is indeed true that the 
thermal diffusion length of the vapor is an order of magnitude greater 
than the diffusion length in the liquid. However, the thermal capacity 
(density times specific heat) of the vapor is two orders of magnitude less 
than that of the liquid. Thus the temperature variations in the liquid 
phase should be an order of magnitude less than the temperature varia- 
tions in the vapor. Therefore, one would expect the principal temperature 
gradients at the interface to lie within the vapor phase rather than the 
liquid. A nonuniform temperature in the vapor would appear to be the 
rule rather than the exception. 

Having raised the specter of a nonuniform, nonequilibrium state inside 
the bubble, one next wonders what to do about it. Clearly this has to be 
considered if improvements are to be made in the theory. Recent attempts 
at this have involved rather gross assumptions and do not appear to be an 
advance on the basic work reported in Professor Plesset’s paper. 

PLESSET (author) : Professor Robertson has asked me to comment on 
the “dissolved gas” effect which is found when a liquid is subjected to 
high pressure. After a liquid such as water has been so treated and when 
it is then tested under ordinary pressure, it shows a large tensile strength. 
This large tensile strength, however, gradually disappears as the liquid 
returns to its normal low tensile strength. I believe that this phenomenon 
is not primarily due to undissolved gas in the usual sense. I think that, 
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rather, it is due to the presence of solid particles which are not wetted by 
the liquid. These particles under ordinary conditions act as holes in the 
liquid and will therefore serve as “nuclei” or sites at which cavities can 
form with moderate tensions. When, on the other hand, the liquid is put 
under high pressure, the liquid will be forced into the surface of the solid 
particles, which will then no longer be unwetted and will no longer act as 
nuclei. Such particles, of course, may have originally stable layers of gas 
adsorbed on their surfaces, but such gas layers would be expected to be 
only a few molecules thick. Upon pressurization, this gas will also be forced 
into the surface of the solid with some liquid and may in this sense be 
dissolved. The gradual return to low tensile strength is to be interpreted 
as the diffusion of liquid out of the solid particles so that it again acts as 
a nucleus. 

While Dr. Hickling has raised some subtle and interesting points in his 
discussion of my paper, I believe that his reservations regarding the 
present treatments of thermal problems in bubble dynamics are not 
justified. He suggests that nonequilibrium behavior in the liquid/vapor 
system will occur even when the interface moves slowly. This question 
has been considered in detail, and it has been shown that the vapor pres- 
sure at a liquid interface will have its equilibrium value (the “Clausius- 
Clapeyron” value) when the motion of the interface is small compared 
with the characteristic vaporization velocity c* where 


(2tt7) 1/2 

Here c is the velocity of sound in the vapor and 7 is the ratio of specific 
heats for the vapor. When the liquid interface is moving into the vapor 
with a velocity R, the pressure in the vapor adjacent to the liquid, p', 
is above the equilibrium value and is given by 


V =P*i 


c*—R 


where p eq is the Clausius-Clapeyron value. When the interface is moving 
away from the vapor, one has 


V' = P", 


c*+R 


These expressions, which are accurate when R is less than c*, show that 
the deviations of p' from p eg are unimportant for interface motions that 
are characterized by small Mach numbers. These slow motions are typical 
of those encountered in most cavitation and boiling situations. 

We can also easily answer the question of the spatial variations in vapor 
pressure throughout a bubble of radius R . The vapor pressure at a bubble 
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wall, as we see from the expressions above, is the equilibrium vapor pres- 
sure corresponding to the liquid temperature at the wall, provided that 
R is small compared to c*. But this condition is a statement that the 
bubble wall has a motion with small Mach number. Such a motion of the 
vapor is essentially subsonic, incompressible motion. In this approxi- 
mation it follows that the pressure in the vapor cannot have any appre- 
ciable spatial variation. The conclusion is that a bubble which has a 
moderate radial velocity, R, has a spatially uniform vapor pressure which 
is given by the equilibrium vapor pressure corresponding to the tempera- 
ture of the liquid interface. 

We cannot yet conclude that the vapor density and vapor temperature 
are uniform throughout the bubble. Although the vapor pressure is 
uniform in the bubble, the vapor density and temperature in the interior 
of the bubble do not necessarily have the values that they assume at the 
boundary. It is clear, however, that both of these latter quantities must 
be spatially uniform if one of them is uniform, since the pressure, density, 
and temperature are connected by the equation of state. We can readily 
see that the vapor temperature is uniform for the moderate values of R 
of present concern. This conclusion follows because the diffusion length 
will then be greater than R, and the temperature in the bubble, while it 
may vary with time, will not have appreciable spatial variation. Dr. 
Hickling is aware that the diffusion length will be large compared with R. 
It is of course well known that the coefficient of thermal diffusion is the 
only physical parameter which appears in the equation for heat flow, and, 
as Dr. Hickling points out, the thermal diffusivity of the vapor is much 
larger than that of the liquid, although the thermal conductivity of the 
vapor is appreciably less than that of the liquid. The only way that the 
thermal conductivity of the vapor or the liquid can enter into a problem, 
however, is through boundary conditions, and the boundary conditions 
should be considered with care since this is where Dr. Hickling is in error 
in his view of the problem. On the liquid side of the interface the heat flow 
per unit time is 


where k L is the thermal conductivity of the liquid, ( dT/dr)i is the tem- 
perature gradient in the liquid at the bubble wall, and 47T.R 2 is the area 
through which the heat flows. The boundary condition on the vapor side 
of the interface is not 

9 

as Dr. Hickling supposes. Rather, the boundary condition in the vapor is 
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L|(frJPp.) 

where L is the latent heat of evaporation or condensation. We see that 
the heat flow into an expanding bubble needed for latent heat of evapora- 
tion is supplied by a temperature gradient in the liquid, but we also see 
that no gradient is developed in the vapor. In the same way, for a shrink- 
ing bubble the latent heat of condensation of vapor is conducted away by 
a suitable gradient in the liquid without the development of a gradient in 
the vapor. Our conclusion is, finally, that a vapor bubble which expands 
or contracts at moderate rates will have a spatially uniform interior, 
uniform in vapor pressure, density, and temperature. 

Dr. Hickling has also made reference to permanent gas bubbles in a 
liquid and again expresses concern that large temperature gradients are 
to be expected in the gas near the bubble boundary . As in the condensible 
vapor case, the thermal diffusivity of permanent gas is much larger than 
that of the liquid, so that the diffusion length may be much larger than 
that of the liquid. The physical situation is straightforward. When the 
diffusion length is large compared with the bubble radius that is, for 
slow processes — the bubble temperature will remain isothermal. The fact 
that the thermal conductivity of the gas is small compared to that of the 
liquid is immaterial, since the temperature gradients in the liquid and in 
the gas are both essentially zero. The boundary condition at the bubble 
wall becomes trivial. The case of rapid motions of a permanent gas bubble 
is perhaps outside the scope of this discussion, but it may have some 
interest nevertheless. If the diffusion length in the gas becomes appre- 
ciably smaller than the bubble radius, the bubble motions will become 
adiabatic in some average sense. So far as the boundary condition in the 
liquid at the bubble wall is concerned, the temperature gradient is 
negligible; no appreciable heat flows into or out of the bubble. Gradients 
develop in the gas, and it is for this reason we must speak of an adiabatic 
behavior in an average sense. 

I hope that these comments represent a clarification of a problem which 
has been beset by difficulties in understanding. 
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